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P-wave heavy-light mesons using NRQCD and D234 
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The masses of S- and P-wave heavy-light mesons are computed in quenched QCD using a classically and 
tadpole-improved action on anisotropic lattices. Of particular interest are the splittings among P-wave states, 
which have not yet been resolved experimentally; even the ordering of these states continues to be discussed in 
the literature. The present work leads to upper bounds for these splittings, and is suggestive, but not conclusive, 
about the ordering. 



1. INTRODUCTION 

The masses of S-wave heavy-light mesons have 
been a valuable testing ground for lattice NRQCD 
studies. The calculations are precise enough to 
test convergence of the NRQCD expansion 
and the effects of quenching|^. Also, experimen- 
tal measurements are precise enough to constrain 
the values of coefficients in the NRQCD action[^. 

For P-waves, the situation is quite different. 
Experiments have not yet resolved the splittings 
among P-wave states, and Isgur[|| has recently 
promoted the suggestion, initiated long ago by 
Schnitzer|Q, that the P-wave states may be in- 
verted relative to the conventional hydrogen-like 
ordering. This inversion is predicted by using 
the experimental K* masses and mixing angle 
as input, then extrapolating from the constituent 
strange quark to charm or bottom via a 1/M- 
expanded nonrelativistic potential. Verification 
of this prediction would thus lend support to 
the treatment of a constituent strange quark as 
"heavy" within a nonrelativistic quark model. 

To date, quenched lattice QCD calculations 
have not seen an inversion |-|^], ahhough the un- 
certainties are often large enough to make firm 
conclusions difficult. 

The present work is an exploration of this 
physics using an improved action on anisotropic 
lattices. The results obtained below agree quali- 
tatively with previous lattice calculations, but a 
quantitative comparison reveals differences. 



2. ACTION 

The lattice action has three terms: gauge ac- 
tion, light quark action and heavy quark ac- 
tion. The entire action is classically and tadpole- 
improved with the tadpole factors, Uo,s and t/o,t, 
defined as the mean links in Landau gauge in a 
spatial and temporal direction, respectively. 

The gauge action includes a sum over 1x2 rect- 
angular plaquettes as well as 1 x 1 elementary pla- 
quettes, and therefore the leading classical errors 
are quartic in lattice spacing. For light quarks, 
a D234 action|H is used with parameters set to 
their tadpole-improved classical values. Its lead- 
ing errors are cubic in lattice spacing. 

The heavy quark action is NRQCD which is 
discretized to give the following Green's function 
propagation 

atHy Ul / atH 
2n J Uo.t V 2ri 

with n = 5 chosen for this work, and the Hamil- 
tonian is truncated at 0{1/M), 
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A tilde denotes removal of the leading discretiza- 
tion errors, leaving the action with quadratic lat- 
tice spacing errors. Notice that the coefficients 
in H are set to their tadpole-improved classical 
values. 
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Table 1 

Heavy-light meson creation operators. 

^^0 (0,/) 
'^1 (0,a,) 

{0,A^<JJ - Aj-aO 

(0,A»Cr» - AjCTj) or 

(0,A»<^j+A,a.O,«^j 

3. RESULTS 

All data presented here come from 2000 gauge 
field configurations on lO'^xSO lattices at /3 = 2.1 
with a bare aspect ratio of as/at = 2, which 
corresponds to at = 0.10 fm. The light quark 
mass is fixed by k = 0.24, which gives rriT^jmp = 
0.52 ± 0.01, ie. rrig ~ mg/l. Fixed time bound- 
aries are used for the light quark propagators so 
they fit naturally into a meson with an NRQCD 
heavy quark propagator. 

Four heavy quark masses have been studied: 
OsM = 1, 3, 5, oo. Calculation of the kinetic mass 
of the heavy-light meson leads to the follow- 
ing charm and bottom quark masses: UgMc ~ 1.2 
and agMb ~ 5. More precise determinations of 
these parameters are not required for the present 
exploratory study. 

Familiar heavy-light meson creation operators 
are used, 

^gt(f)f](f)r(f)q(x), (4) 

X 

where fl{x) is given in table |l| and the smearing 
operator is 

r(f) = [l + c,A(2)(f)]"=. (5) 

All plots shown here use {cs,ns) = (0.15,10) at 
the source (which is fixed to timestep 4) and a 
local sink. 

Fig. |l] shows the ^5*0 and -^Pi simulation en- 
ergies in the bottom region. In NRQCD, only 
energy differences are physical but this plot is an 
indication of the data quality. As for all plots in 
this work, the uncertainties are determined from 
5000 bootstrap ensembles. It should be noted 
that the ^Pi-^Sq splitting is in agreement with 
previous lattice determinations (eg. 



Fig. I shows the ^Si-^So mass splitting and 
the kinetic shift of the ^Sq simulation energy in 
the bottom region. In physical units, the ^Si-^Sq 
splitting is found to be 34±5 MeV, in agreement 
with the known quenched result (eg. 0J^-0]). 

Fig. H displays the mass splittings among P- 
wave states at all available quark masses. The 
•^Pi is not shown since it is not distinguishable 
from the ^Pi. Near the source, the states are 
organized according to the familiar hydrogen-like 
pattern: Pq, Pi, P2 from lightest to heaviest. 

As time increases, the uncertainties grow such 
that one might be tempted to define a plateau 
which begins rather close to the source. However, 
the central values of the mass splittings tend to 
decrease for increasing time, so a negative mass 
splitting (i.e. an inverted spectrum) cannot be 
ruled out conclusively. 

Using at = 0.10 fm. Fig. || indicates that 
the '^P2-'^Po splitting for bottom mesons is less 
than 100 MeV, which should be contrasted with 
Ref. §, where lattice NRQCD gave 183±34 MeV. 
While it is true that the two determinations stem 
from different methods (different lattice actions, 
anisotropic versus isotropic lattices, different light 
quark masses, slightly different temporal lattice 
spacings, ...), it is disconcerting that the final 
results are not in better agreement. 

In conclusion, this work has led to an upper 
bound for the '^Pj-'^Po (bottom) splitting which 
lies below a previous lattice determination. An 
inverted spectrum, although not suggested by the 
calculation, cannot yet be ruled out. 

a M = 5 
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Figure 1 . The ^ 5*0 and ^ Pi simulation energies in 
lattice units at QsM — 5. 
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Figure 2. The S-wave mass splittings in lattice 
units at OsM = 5. The nonzero momentum "p" 
is 27r/(10as). 



E 0.0 




0.2 



E 0.0 



ml 



15 



E 0.0 



It 



15 



static heavy quark 



E 0.0 - 




Figure 3. The P-wavc mass splittings in lattice 
units. 



